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We study the interface between a fractional topological insulator and an ordinary insulator, both 
described using holography. By turning on a chemical potential we induce a finite density of matter 
localized at the interface. These are gapless surface excitations which are expected to have a 
fermionic character. We study the thermodynamics of the system, finding a symmetry preserving 
compressible state at low temperatures, whose excitations exhibit hyperscaling violation. These 
results are consistent with the expectation of gapless fermionic excitations forming a Fermi surface 
at finite density. 



I. INTRODUCTION 

Holography is a fairly recent approach to the study 
of strongly interacting systems which have gravitational 
duals. As such, it is ideal to model and gain qualitative 
understanding of otherwise ill-understood phenomena in 
condensed matter systems. The approach has been ap- 
plied successfully to the study of homogeneous phases of 
holographic matter, and has provided insight into diverse 
phenomena such as superconductivity, quantum phase 
transitions, and non-Fermi liquid behavior. 

Many interesting phenomena in modern condensed 
matter physics involve boundary excitations at the in- 
terface of two different bulk materials. Under suitable 
circumstances, those boundary excitations are gapless 
by virtue of topological considerations. Those consid- 
erations protect the excitations from developing a gap 
in the presence of noise and other small perturbations. 
These localized excitations are particularly interesting 
when they are predicted to be purely fermionic, as they 
are a natural candidate for an holographic non-Fermi liq- 
uid p. 

In this note we initiate the study of gapless surface 
excitations living on interfaces of different holographic 
materials. To that end, we consider the holographic 
model constructed in [2] (for related development see [3J- 
8 ) which describes the interface of a fractional topologi- 
cal insulator and an ordinary (topologically trivial) insu- 
lator. Briefly, the model describes a fermion coupled to 
massless gauge fields, with spatially varying mass profile 
m{x); index theorems then guarantee the existence of a 
fermionic bound state on the surface 1 . Since the fermions 
are expected to be strongly correlated (by virtue of their 
coupling to deconfined gauge fields), it is interesting to 
study the physics of these boundary excitations using 
various probes. 



The model also has charged bosonic excitations, but they are not 
expected to form bound states at the interface. We find below 
that there is no symmetry breaking in the model, i.e. the bosons 
are not condensed. 



Here we start such an investigation by studying the 
thermodynamics of the system. We find that the matter 
localized at the interface is compressible - the density 
d varies smoothly with the chemical potential \x, with 
d'(p) =^ for all values of the chemical potential (and 
T <C (J.). Furthermore, the excitations around this state 
seem to exhibit a scaling symmetry with hyperscaling 
violation, within the precision of our numerics. These 
results suggest the existence of a Fermi surface for the 
boundary excitations. 

The plan of this note is as follows: in section II we 
introduce the holographic setup and describe the inho- 
mogeneous backgrounds we construct, where the surface 
excitations are put at finite density. In section III we 
discuss the thermodynamics of the system. We end by 
discussing directions for further exploration of the nature 
of these and other boundary excitations. 

II. THE HOLOGRAPHIC SETUP 

We are interested in describing the surface excitations 
at non-zero temperature T and finite chemical potential 
p. The starting point is then a black hole in AdS§ x S 5 . 
We use the conventions of [SHU], work in the Poincare 
coordinates, and set the radius of curvature of AdS$ to 
unity: 

ds 2 = \ 2 T 2 (~Q^dt 2 + h{p)dxA +^- + dfll 
2 V h{p) J p 2 

dQj = d6 2 + sin 2 9 dflj + cos 2 9 dej) 2 
fiP)- 1 -^ h{p) = l + ^ 

We are interested in 7-branes embedded in AdS§ and 
wrapping the 3-sphere inside S 5 , localized at an angle 9 
on the five-sphere. The 7-brane embedding is character- 
ized by the variation of 9 as function of a radial coor- 
dinate p in AdS, and in our case one of the boundary 
coordinates x. The embedding is a minimal surface sub- 
ject to asymptotic boundary conditions which encode the 
parameters of the problem: temperature, chemical poten- 
tial and the mass parameter m. The precise definition of 
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those asymptotic quantities appears below. 

The embedding can end smoothly when the 3-sphere 
shrinks to zero size [T^], psinO — 0, which can happen in 
one of two ways: 

• If p 7^ and s'm9 — 0, the 7-brane ends at 
a finite radial coordinate in AdS, away from the 
horizon, whose value depends (monotonically) on 
the mass parameter m. Such embeddings are 
dubbed " Minkowski embeddings" ; they correspond 
to states with no density 2 . For zero temperature 
and p < \m\ these are the only possible embeddings 
and therefore in this parameter range the system is 
necessarily at zero density. 

• On the other hand, if p = 0, then —1 < sinf? < 1 
parametrizes the inner boundary of the probe brane 
which lies on the black hole horizon. Such embed- 
dings, called "black hole" embeddings [TUHH], can 
support finite density configurations and are possi- 
ble, at zero temperature, only when p > \m\. 

In order to describe an interface we study the system 
with a mass profile m(x) depending on a single boundary 
direction x, interpolating between m{x) — —M, a case 
which by convention we call the topological insulator, and 
m{x) = M which corresponds in our sign conventions to 
a regular insulator. The system is insulating in the bulk, 
but supports gapless excitations on the interface (where 
m(x) = 0). We are interested in probing the physics of 
those excitations. 

To that effect we turn on a constant 3 chemical po- 
tential p. The effect of the chemical potential is easiest 
to understand at zero temperature. In that case, if we 
choose to satisfy p < M, the chemical potential would in- 
duce no finite density in either one of the bulk insulators, 
but it would induce finite density at the interface, in the 
region where m{x) < p. The embedding then interpo- 
lates between the Minkowski embedding (ending at some 
finite radial coordinate), and the black hole embedding 
(ending at the black hole horizon). 

In practice, in order to avoid numerical instabilities 
associated with the abrupt transition between the two 
types of embeddings 4 , we choose the mass profile such 
that a small density is induced at the edges of the spa- 
tial interval. While the embedding is then of the " black 
hole" type throughout the spatial interval, the quantities 
calculated arc largely independent of the mass profile, as 
long as the density induced is at least moderately peaked 
at the interface. 



2 Minkowski embeddings with strings attached are possible in prin- 
ciple, but they give subleading contributions to the thermody- 
namics for any finite density. 

3 For a system in diffusive thermal equilibrium the chemical po- 
tential is constant even for inhomogeneous configurations. 

4 Though it is not necessary for our purposes, it is interesting to 
construct solutions which correspond to phase boundaries in the 
holographic context, along the lines of |13| . 



Ansatz and Equations of Motion 

To describe the embedding we choose to parametrize 
the worldvolume in terms of the radial coordinate p and 
the spatial coordinate x. In our parametrization, the em- 
bedding is characterized by a single function x(p, x) = 
cosf?(p, x). The asymptotic behavior of the coordinate 
X at large p, x(P: x ) ~~ ^ , determines the mass pro- 
file 771(2;), we choose to approximate a step function as 
m(x) — 1+ 2 e M ax — M for a> 1. As long as the mass pro- 
file is sufficiently narrow and steep, the parameters a, M 
do not significantly influence the physical quantities we 
calculate 5 . 

Having chosen a parametrization, we add a worldvol- 
ume gauge field Aq (choosing 2ira' = 1). It is straight- 
forward to find the following DBI action: 

Sdbi = 7VT 4 p 3 (1 - X 2 ) / h y/I x + I Aq + I int 

J* = i-x 2 + (3 P x) 2 + ^f 

2 ( 2h{d p A,f 2{d x A f \ 

u - { 1 - x ) p pi p J 

lint = -p 2 (d p xd x A - d p A d x x) 2 

where the form of I x and Ia reveals the elliptic nature of 
the equations. The field Aq and the coordinate x are cho- 
sen to be dimensionless, like all other quantities appear- 
ing in the action. With this choice we measure dimen- 
sionful quantities in units of the temperature T. There- 
fore, the asymptotic value of Aq at large p determines 
the chemical potential in units of T, Aq — > p = ^ . 

The regularity of the solution near the boundary of the 
embedding is discussed in [T3]. For a black hole embed- 
ding we need to satisfy d P x = and Aq = at p = 0. 
Asymptotically in the p direction, we set a UV cutoff 
p = A and choose A (p = A, x) = p to set the chem- 
ical potential, and x(P — A, x) — to specify the 
mass profile. Finally, we choose homogeneous Neumann 
boundary conditions for both functions at the boundaries 
of the spatial interval. 

The Solutions 

We solve the equations of motion numerically using 
pseudospectral methods implemented in Matlab, utiliz- 
ing a Chebyshev grid to discretize the equations, and 
solving the resulting set of algebraic equations using 
a Newton-Raphson-Kantorovich iteration. The conver- 
gence of the solution with the linear size of the grid is 
exponential. The results shown are for a grid of 900 



J We checked that changing those parameters by two orders of 
magnitude changes the calculated quantities by less than one 
part in 10 4 . 
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sity. 
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FIG. 1: The embedding function \ near the horizon. 
The physical parameter is (m/T = 2. 



points (giving accuracy of about one part in 10 6 for the 
values of the fields). 

Before extracting physical properties, we briefly dis- 
cuss the features of the solutions. In figures (1-3) we 
display a few characteristic solutions, at various values 
of the parameters. 

In figure 1 we display the embedding function x- The 
spatial variation of x increases towards the horizon, 
where its profile is correlated with the boundary den- 
sity: the density is maximal for x — and decreases to 
zero as x ~ ^ il- The value % = ±1 at the horizon would 
correspond to the transition to the Minkowski embedding 
which carries no density. 

In figure 2 we show a contour plot of the bulk gauge 
field. The spatial variation is small around the horizon 
and grows towards the boundary. Note that the leading 
behavior of the gauge field near the boundary is spatially 
homogeneous, and any inhomogeneity is manifested in 
the behavior of the subleading term, related to the den- 





FIG. 2: Contour plot of the gauge potential Aq that 
shows the spatial inhomogeneity. The horizontal 
direction is a section of the the radial coordinate. The 
physical parameter is /j/T = 3.4. 



III. THERMODYNAMICS 

The thermodynamics of the system is completely spec- 
ified when expressing the number density n in terms of 
the chemical potential \x and the temperature T. Equiv- 
alently one can specify the free energy F(/j,, T), then 
ri = We choose to calculate the density n, which 

is read off the subleading fall-off of the gauge field, 
Aq — > jl — -^2- n(fl) as p — > oo, where ft, = ^, 

The subleading term in the expansion of Xi related to 
the chiral condensate in e.g. [10], is found to vanish at the 
interface for all our solutions. This is a result of a parity 
symmetry preserved by our mass profile, and indicates 
the absence of symmetry breaking and the existence of 
non-trivial fluid at the interface. The physics of that fluid 
depends only on mu and T, leading to the dimensional 
analysis below. 

Since our background is inhomogeneous, there is sep- 
arate information in densities per unit volume (at fixed 
spatial position x), and the corresponding densities per 
unit area (when integrating over the spatial direction x). 
Here we mainly discuss the free energy and number den- 
sities per unit volume, at the center of the interface x — 0, 
ending with comments on the integrated quantities. 

Based on dimensional analysis in 3+1 dimensions, F = 
T 4 J-(jl), since T, \x are the only scales in the problem [15] . 
For relativistic systems which form a Fermi surface at low 
temperatures, for example for free massless fermions 6 , 
T(Ji) ~ /i 4 + .... The leading term in this expansion 
corresponds to n ~ /i 3 , the contribution to the density 
from a putative Fermi surface. 

The subsequent terms in the small /j expansion give 
information on low energy excitations of the system. In 
particular, the second term in the expansion gives the 
leading temperature dependence of the free energy at low 
temperatures. Thus, this term is responsible for the tem- 
perature scaling of the entropy S in that regime. The 
entropy scales differently for systems with or without a 
Fermi surface. For degenerate relativistic fermions, one 
has Sf^TaX low temperatures, whereas for relativistic 
bosons Sb ~ T 3 (see for example [E]). The larger en- 
tropy at low temperatures reflect the fact that systems 
with a Fermi surface have more low energy excitations in 
the degenerate limit 7 . 

The difference between bosons and fermions can also 
be parametrized in terms of the hyperscaling violation 



6 The normalization of this expression depends on the normaliza- 
tion of the fields in the microscopic theory, for example for us 
this would involve the normalization of the fermion bound state 
wavefunction. 

7 Note that the scaling of the entropy depends only the density 
of low energy states and does not assume those excitations are 
weakly interacting. 
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FIG. 3: Interface density (normalized to unity at p, = 1) 
as function of the chemical potential, on a log-log scale. 
The power law behavior of n{Ji) at large jl is apparent, 
as are deviations from this behavior at lower values of jl. 



confidence 9 , as c = (the value required for hyperscaling 
to hold) is about 3 sigma away from the central value 
for the fit for c. This fits the expectation of hyperscaling 
violation exponent 9 — 2, the value corresponding to 
small excitations of a Fermi surface. 

In summary, we have examined the density per unit 
volume induced on the interface, revealing structure con- 
sistent with the fermionic nature of the boundary excita- 
tions. Unfortunately, the densities per unit area cannot 
be reliably computed using our solutions. Tracking the 
change in the boundary layer width as function of the pa- 
rameter jl is necessary for that purpose. However, in our 
appraoch the width of the domanin wall is determined by 
a rather than by dynamical properties of the interface. 
It would be interesting to improve on our discussion and 
calculate reliably the densities per unit area; we leave 
that for future work. 



exponent 9, such that S ~ T 3 ~ e with 9 = for bosons 
and 9 = 2 for fermions. The origin of the difference 
in scaling is the nature of the excitations of the Fermi 
surface. The low energy excitations are constrained to 
be close to the Fermi surface, thus reducing the effective 
dimensionality of the system. 

After parametrizing the possible small temperature be- 
havior, we can present the results of our model. Figure 
(4) displays a log-log plot of the density as function of the 
chemical potential (both measured in units of tempera- 
ture) . The small temperature expansion discussed in this 
section corresponds to the large jl limit. The functional 
dependence encoded in figure (4) reveals a few interesting 
features of the system studied. 

The most robust feature is the leading behavior at large 
jl, n(jl) ~ ji a , with a ~ 3 to numerical accuracy. De- 
pending on fit options, both for reading off the density 
for each solutions and for fitting the curve n(ji), a is fairly 
robust, changing between 2.95 and 3.1 (for the displayed 
graph a = 2.98). 

This result is consistent with the expectation that the 
ground state of the system is a Fermi surface. This re- 
lation implies that the derivative of the charge density 
d with respect to the chemical potential is non-zero at 
low temperatures, i.e. the state at the interface is com- 
pressible. A Low temperature compressible state, in the 
absence of symmetry breaking, can be taken as indicative 
of a Fermi surface (see e.g. the discussion in [17 ). 

Reading off the subleading terms in the jl expansion is 
less robust. Nevertheless, one can develop an asymptotic 
expansion in 4, a low temperature expansion. To probe 
nature of the low energy excitations, parametrized by the 
value of 9, we fit the asymptotic data to the form 8 n(/2) ~ 
ji +cjl. We find that the value for c is non-zero with great 



IV. CONCLUSIONS AND OUTLOOK 

In this note we constructed inhomogeneous solutions 
describing the holographic dual of an interface between a 
fractional topological insulator and an ordinary insulator. 
We have used these solutions to start the investigation of 
the gapless strongly correlated fermions expected to live 
on that interface. Indeed, by studying the low tempera- 
ture expansion of the thermodynamics of the system we 
found results consistent with these expectations. 

A few directions for further study are immediately ap- 
parent. It is interesting to study more general models 
with fermions bound to an interface, for example mod- 
els with a non-trivial dynamical exponent z, or models 
with fermionic excitations only - though bosons are not 
expected to form gapless bound states at the interface, 
they are likely to limit the precision of the analysis. It is 
also desirable to study such systems at zero temperature 
and finite density. In that limit the effects of a degen- 
erate Fermi gas are likely to be more pronounced and 
can be more cleanly studied. As mentioned in the text, 
moving away from the "thick" domain wall regime is an 
interesting, albeit challenging, problem. 

Most interestingly, since the boundary fermions are ex- 
pected to be strongly correlated (i.e. coupled to massless 
gauge fields), these types of systems are natural candi- 
dates for an holographic dual for a non-Fermi liquid. It is 
therefore particularly interesting to search for further sig- 
natures of the fermionic nature of the interface, for exam- 
ple oscillations in response functions, or non-analyticity 
in the fermionic spectral densities of the system. We ex- 
pect that general enough models will exhibit characteris- 
tics of strange metallic behavior, and we hope to return 



A quadratic term would correspond to 9 = 3, violtaing the bound 
9 < d — 1, based on the asymptotic scaling of the entanglement 
entropy expected to hold in all local quantum field theories |18| . 



Study of subsequent terms in the jl expansion should reveal con- 
tributions from bosonic degrees of freedom. However, the sce- 
nario of having only bosonic contributions is excluded. 
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to some of these questions in the near future. 

Finally, beyond that particular example of an holo- 
graphic interface, there is a whole range of fascinating 
interface physics to be explored holographically. We ex- 
pect that the methods used here will be useful in explor- 
ing this research direction. 



eel Franz, Jeremy Heyl, Gordon Semenoff, Darren Smyth, 
Jared Stang and Evgeny Sorkin and especially from a cor- 
respondence with Kristan Jensen and Andreas Karch, for 
which I am grateful. 



Acknowledgments 



I am supported by a discovery grant from NSERC of 
Canada. I have benefitted from conversations with Mar- 



[1] See comments in Andreas Karch's KITP blackboard talk, 
http://online.kitp.ucsb.edu/online/ adscmt 11/ karch / 

[2] C. Hoyos-Badajoz, K. Jensen and A. Karch, "A Holo- 
graphic Fractional Topological Insulator," Phys. Rev. D 
82, 086001 (2010) [arXiv:1007.3253| [hep-th]]. 

[3] A. Karch, "Electric-Magnetic Duality and Topologi- 
cal Insulators," Phys. Rev. Lett. 103, 171601 (2009) 
|arXiv:0907.l528| [cond-mat.mes-hall]]. 

[4] J. Maciejko, X. -L. Qi, A. Karch and S. -C. Zhang, "Frac- 
tional topological insulators in three dimensions," Phys. 
Rev. Lett. 105, 246809 (2010) [arXiv:1004.3628| [cond- 
mat.str-el]]. 

[5] A. Karch, J. Maciejko and T. Takayanagi, "Holographic 
fractional topological insulators in 2+1 and 1+1 dimen- 
sions," Phys. Rev. D 82, 126003 (2010) [arXiv:1009.299T| 
[hep-th]]. 

[6] S. Sun and A. Karch, "Finite size effects on the statisti- 
cal angle of an electron induced dyon in proximity to a 
topological insula tor," Phys. Rev. B 84, 195115 (2011) 
[arXiv:1108."0700| [cond-mat.mes-hall]]. 

[7] J. Maciejko, X. -L. Qi, A. Karch and S. -C. Zhang, 
"Models of three-dimensional fractional topological in- 
sulators," arXiv:1111.6816 [cond-mat.str-el]. 

[8] B. Swingle, M. Barkeshli, J. McGreevy and T. Senthil, 
"Correlated Topological Insulators and the Fractional 
Magnetoelectric Effect," Phys. Rev. B 83, 195139 (2011) 
[arXiv: 1005.10761 [cond-mat . str-el] ] . 

[9] D. Mateos, R. C. Myers and R. M. Thomson, -'Holo- 
graphic phase transitions with fundamental matter," 
Phys. Rev. Lett. 97, 091601 (2006) |hep-th/0605046 . 



[10] S. Kobayashi, D. Mateos, S. Matsuura, R. C. Myers 
and R. M. Thomson, "Holographic phase transitions at 
finite baryo n density," JHEP 0702, 016 (2007) hep- 
th/0611099] . 

[11] D. Mateos, S. Matsuura, R. C. Myers and R. M. Thom- 
son, "Holographic phase transitions at finite chemical po- 
tential," JHEP 0711, 085 (2007) |arXiv:0709.1225l [hep- 
th]]. 

[12] A. Karch and E. Katz, "Adding flavor to AdS / CFT," 
JHEP 0206, 043 (2002) hep-th/0205236 . 

[13] O. Aharony, S. Minwalla and T. Wiseman, "Plasma- 
balls in large N gauge theories and localized black holes," 
Class. Quant. Grav. 23, 2171 (2006) |hep-th/05072T9] . 

[14] V. P. Frolov, "Merger Transitions in Brane-Black-Hole 
Systems: Criticality, Scaling, and Self-Similarity," Phys. 
Rev. D 74, 044006 (2006) [gr-qc/0604114] . 

[15] F. Nogueira and J. B. Stang, "Density versus chemical 
potential in holographic field theories," Phys. Rev. D 86, 
026001 (2012) [arXiv: 1111. 2806] [hep-th] ] . 

[16] L. Homorodean, "Thermal Properties of Degenerate Rel- 
ativistic Quantum Gases." Modern Physics Letters B, 
Volume 17, Issue 18, pp. 991-1010 (2003). 

[17] L. Huijse and S. Sachdev, "Fermi surfaces and gauge- 
gravity duality," Phys. Rev. D 84, 026001 (2011) 
[arXiv:1104.5022| [hep-th]]. 

[18] L. Huijse, S. Sachdev and B. Swingle, "Hidden Fermi 
surfaces in compressible states of gauge-gravity duality," 
Phys. Rev. B 85, 035121 (2012) |arXiv:1112.0573"l [cond- 
mat. str-el]]. 



